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Abstract. We discuss the nature of phase transitions in the self-gravitating Fermi gas 
at non-zero temperature. This study can be relevant for massive neutrinos in Dark 
Matter models and for collisionless self-gravitating systems experiencing a "violent 
relaxation" . Below a critical energy (in the microcanonical ensemble) or below a critical 
temperature (in the canonical ensemble) , the system undergoes a gravitational collapse 
leading to a compact object with a massive degenerate nucleus (fermion ball). 



1 Introduction 

The self-gravitating Fermi gas model was introduced by Fowler in 1926 in order 
to explain the puzzling nature of white dwarf stars. When a star has exhausted 
its nuclear fuel, it undergoes a gravitational collapse until gravity is balanced 
by the pressure of a degenerate gas of electrons (on account of Pauli's exclusion 
principle). This model was later on improved by Chandrasekhar (1930) who took 
into account special relativity effects in the degenerate equation of state of elec- 
trons and discovered that white dwarf stars have a limiting mass of ~ 1.4 M Q . 
Then, in 1939, Oppenheimer & Volkoff considered the extension of this model to 
the framework of general relativity in connexion with the structure of neutron 
stars. They found again a limiting mass of ~ 0.7 M Q above which the pres- 
sure of the neutrons cannot sustain gravity anymore. In that case, the star is 
expected to collapse into a black hole. In these studies, the Fermi gas is com- 
pletely degenerate since the thermal energy kT is much smaller than the Fermi 
energy. The self-gravitating Fermi gas at finite temperature was investigated by 
Hcrtel & Thirring (1971). They showed that below a critical temperature a ho- 
mogeneous gas of fermions undergoes a first order phase transition leading to 
a compact object with a degenerate core. Application of these results to Dark 
Matter was considered by Bilic & Viollier (1997). They proposed that Dark 
Matter could be made of massive neutrinos (with m ~ 15keV) in equilibrium 
with a radiation background imposing its temperature. In the langage of statis- 
tical mechanics, this corresponds to the canonical ensemble. By cooling below a 
critical temperature, a condensed phase emerges consisting of a quasidegenerate 
"fermion ball". These fermion balls may provide an alternative to black holes 
that are reported to exist at the center of galaxies. At large distances, the r~ 2 
law of density decrease of an isothermal gas is consistent with the flat rotation 
curves of spiral galaxies. This could provide an attractive model of Dark Matter 
but the dynamical mechanism leading to a "fermion ball" remains to be clearly 
identified. 
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The self-gravitating Fermi gas model also appeared in a completely different 
context, indcpcndant of quantum mechanics. In 1967, Lynden-Bcll argued that 
a self-gravitating system far from mechanical equilibrium would rapidly relax 
towards a virialized state, due to the strong fluctuations of the gravitational 
potential. Since this process of violent relaxation is essentially collisionless, the 
coarse-grained distribution function / cannot increase and it must satisfy an 
effective exclusion principle / < 770 , where 770 is the maximum value of the initial 
distribution function. This upper bound is a consequence of the Liouville theo- 
rem. Assuming ergodicity, Lynden-Bcll predicted that / should converge towards 
a Fermi-Dirac distribution, or a superposition of Fermi-Dirac distributions. Re- 
laxation equations towards these maximum entropy states (on a coarse-grained 
scale) were proposed by Chavanis, Sommeria and Robert (1996). In that con- 
text, the proper thermodynamical ensemble is the microcanonical ensemble since 
energy is conserved during the course of the evolution. Fermi-Dirac spheres were 
computed by Chavanis & Sommeria (1997). They found a wide variety of nu- 
clear concentration depending on the degree of degeneracy and on the value of 
energy. For high energies, the system is almost uniform. For intermediate ener- 
gies, the core is partially degenerate. For low energies, the maximum entropy 
(i.e. most probable) state consists of a massive degenerate nucleus surrounded 
by an isothermal halo. They argued that degeneracy (in Lynden-Bell's sense) 
could stabilize the system and play a role in galactic nuclei and Dark Matter. 
It is however unclear whether violent relaxation can lead to massive "fermion 
balls" . Indeed, it is in general advocated that the fluctuations of the gravita- 
tional potential fade before the system has developed high density contrasts. To 
our point of view, this remains a matter for further investigation. 

Recently, we have carried on a detailed study of phase transitions in the self- 
gravitating Fermi gas (Chavanis 2002). We showed that this system exhibits a 
rich structure with the occurence of three types of phase transitions of zeroth, 
first and second order. We worked at a general level without specifying the 
source of degeneracy (quantum mechanics or violent relaxation) and we described 
the complete structure of the equilibrium phase diagram, for arbitrary values 
of control parameters and arbitrary degree of degeneracy. We worked both in 
the microcanonical and canonical ensembles, emphasizing the inequivalence of 
these ensembles for long-range systems. The description of the microcanonical 
ensemble and the relation between the structure of Fermi-Dirac spheres and 
classical isothermal spheres investigated by Antonov (1962) and Lynden-Bell & 
Wood (1968) is a specificity of our approach. In the following, we give a short 
summary of this study. 



2 Phase transitions in self-gravitating systems 

We consider a system of N fermions of mass m interacting via Newtonian gravity. 
At statistical equilibrium, the system is described by the Fermi-Dirac distribu- 
tion 

/= 7^—7, (P = m/kT) (1) 

1 + Xe^+^ 
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coupled to the Poisson equation 



A$ = AtiG I fd 3 v 



(2) 



This determines a self-consistent meanfield equation for the gravitational poten- 
tial <P. For spherically symmetrical systems, this is just an ordinary differential 
equation, which can be solved numerically by usual means. Then, the fugac- 
ity A -1 > and the inverse temperature (3 can be related to the total mass 
M and total energy E of the system. In Eq. (Q), 770 represents the maximum 
allowable value of the distribution function, i.e. / < ryo- For a quantum gas, 
r/o = (2s + l)m 4 /h 3 where s is the spin of the particles. In the fully degenerate 
limit / ~ rjo, the system is equivalent to a polytrope of index n = 3/2. In the 



The Boltzmann-Poisson system (|3|) has been studied in relation with the struc- 
ture of isothermal stellar cores and globular clusters. It is well-known that the 
density of an isothermal gas decreases at large distances like r~ 2 . Therefore, the 
total mass of the configuration is infinite and we need to introduce truncated 
models (Chavanis 1998) or confine the system within a box of radius R. 




(3) 
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Fig. 1. Equilibrium phase diagram for classical isothermal spheres. For A > A c or 
77 > ?i c , there is no hydrostatic equilibrium and the system undergoes a gravitational 
collapse. 
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Fig. 2. Equilibrium phase diagram for self-gravitating fermions with different values 
of degeneracy parameter p. For p — > +00, the spiral makes several rotations before 
unwinding pH. 



In the non degenerate limit, the equilibrium phase diagram (E,T) is rep- 
resented in Fig. [I]. The curve has a striking spiral behaviour parametrized by 
the density contrast 1Z = p(0)/p(R) going from 1 (homogeneous system) to +00 
(singular sphere) as we proceed along the spiral. There is no equilibrium state 
below E c = -0.335GA/ 2 /i? or T c = ^MeR • In that case > the s y stem is expected 
to collapse indefinitely. This is called gravothermal catastrophe in the micro- 
canonical ensemble (fixed E) and isothermal collapse in the canonical ensemble 
(fixed T) . Dynamical models show that the collapse is self-similar and develops a 
finite time singularity (see Chavanis et al. 2002 and references therein). However, 
although the central density goes to +00, the shrinking of the core is so rapid 
that the core mass goes to zero. Therefore, the singularity contains no mass and 
this process alone cannot lead to a black hole. 

It is also important to stress that the statistical ensembles are not inter- 
changeable for systems with long-range interaction, like gravity. In the micro- 
canonical ensemble, the series of equilibria becomes unstable after the first turn- 
ing point of energy (MCE) corresponding to a density contrast of 709. At that 
point, the solutions pass from local entropy maxima to saddle points. In the 
canonical ensemble, the series of equilibria becomes unstable after the first turn- 
ing point of temperature (CE) corresponding to a density contrast of 32.1. At 
that point, the solutions pass from minima of free energy (F = E — TS) to 
saddle points. It can be noted that the region of negative specific heats between 
(CE) and (MCE) is stable in the microcanonical ensemble but unstable in the 
canonical ensemble as expected on general physical grounds. 
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3. Equilibrium phase diagram for Fermi-Dirac spheres with a degeneracy param- 
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(GEM) for A < At(^) and local entropy maxima (LEM), i.e. metastable states, for 
A > At([i). Points C form the "condensed" phase. They are LEM for A < A t (n) and 
GEM for A > A t (fi). Points B are unstable saddle points (SP) and contain a "germ". 
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Fig. 4. Entropy of each phase versus energy for fj, — 10 s . A first order phase transition 
is expected at At{n) at which the two stable branches (solutions A and C) intersect. 
However, the entropic barrier B probably prevents this phase transition |8|]l4|. 
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When degeneracy is taken into account, the structure of the equilibrium phase 
diagram depends on the value of the degeneracy parameter [i = i] V512nHPMB^ 
The classical limit is recovered for /i — * +00. We see that the inclusion of de- 
generacy has the effect of unwinding the spiral (Fig. ||). For large values of 
the degeneracy parameter, the equilibrium phase diagram is depicted in Fig. 
[| The solutions on the upper branch (points A) are non degenerate and have 
a smooth density profile; they form the "gaseous phase" . The solutions on the 
lower branch (points C) have a "core-halo" structure with a massive degenerate 
nucleus and a dilute atmosphere; they form the "condensed phase". The den- 
sity profiles of these solutions are given in Ref . || . By comparing their entropy 
(Fig. |]), we would expect that a first order phase transition from the gaseous 
phase to the condensed phase occurs at the transition energy Et(n). This is, 
however, unlikely because the probability of a fluctuation able to induce this 
phase transition is extremely weak Therefore, the metastable gaseous states 
with E < Et are probably physical for the time scales involved in astrophysical 
situations. In any case, a phase transition must occur at the critical energy E c at 
which the gaseous branch disappears. Below that energy, the system undergoes 
a gravothermal catastrophe but, for self-gravitating fermions, the core ceases 
to shrink when it becomes degenerate. Since this collapse is accompanied by 
a discontinuous jump of entropy, this is sometimes called a zeroth order phase 
transition. The resulting equilibrium state (point D) possesses a small degen- 
erate nucleus which contains a moderate fraction of the total mass (~ 20% for 
\i = 10 5 ). The rest of the mass is diluted in a hot envelope held by the box. In 
an open system, it would be dispersed at infinity. 
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Fig. 5. Equilibrium phase diagram for Fermi-Dirac spheres with /i = 10 3 . 
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Fig. 6. Enlargement of the phase diagram near the tricritical point in the canonical 
ensemble. The Maxwell construction determining the transition temperature r]t(^) is 
made explicitly (dashed areas). For ficTP = 82.5 the Maxwell plateau disappears and 
the A — r\ curve presents an inflexion point at Actp — 0.5, rjcrp — 3.06. At that point, 
the specific heat becomes infinite and the transition is second order. This diagram is 
remarkably similar to the liquid/gas transition for an ordinary fluid. 

For smaller values of the degeneracy parameter, the equilibrium phase di- 
agram is represented in Fig. ^[ The curve rj(A) is now univalued and the first 
order phase transition in the microcanonical ensemble is suppressed: all the equi- 
librium states are global entropy maxima. For large energies, they are almost 
homogeneous and for smaller energies they have a "core-halo" structure with a 
partially degenerate nucleus. As energy is progressively decreased, the nucleus 
becomes more and more degenerate and contains more and more mass until a 
minimum energy E min = -2.36 G 2 M 7 / 3 77o /3 at which all the mass is in the 
completely degenerate nucleus of radius i?* = 0.181 ?/ 2 ' 3 G _1 M -1 / 3 . In that 
case, the system has the same structure as a cold white dwarf star. Therefore, 
depending on the degree of degeneracy and on the value of energy, a wide variety 
of nuclear concentrations can be obtained in the Fermi gas. 

We can also consider the situation in which the system is in contact with a 
heat bath which imposes its temperature (canonical description). Considering 
again the case fj, = 10 3 in Fig. ||, we note that the curve A{rj) is multivalued. 
We expect therefore a first order phase transition to occur at a transition tem- 
perature T t determined by a Maxwell construction. This phase transition should 
be accompanied by a huge release of latent heat. This may not be physically 
realizable and the true collapse will rather occur at T c . The outcome of this 
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Fig. 7. Same as Fig. |6| near the tricritical point in the microcanonical ensemble 
(fi MT P ~ 2600, A MTP ~ 0.38, vmtp ^ 1-68). 



collapse is the formation of a fermion ball containing almost all the mass. Con- 
sidering Fig. § again, we see that the phase transitions are suppressed in the 
microcanonical ensemble for fj, < [amtp — 2600 and in the canonical ensem- 
ble for /j, < hctp — 82.5. At these tricritical points, the two phases (gaseous 
and condensed) merge. This characterizes a second order phase transition (see 
Figs. |6|-(7]). The analogies and the differences with the liquid/gas transition are 
described in Ref. |J. 
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